Model for the free- volume distributions of equilibrium fluids 
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We introduce and test via molecular simulation a simple model for predicting the manner in which interpar- 
ticle interactions and thermodynamic conditions impact the single-particle free-volume distributions of equilib- 
rium fluids. The model suggests a scaling relationship for the density-dependent behavior of the hard-sphere 
system. It also predicts how the second virial coefficients of fluids with short-range attractions affect their 
free-volume distributions. 



I. INTRODUCTION 

Liquid-state theory aims to provide a framework that links 
the interparticle interactions of a fluid with its local structure, 
thermodynamic properties, and transport coefficients. One of 
the central quantities for characterizing the structural order 
is the static structure factor S{k) (or, equivalently, the pair 
correlation function g{r))} The structure factor can be read- 
ily measured by scattering experiments, computed via molec- 
ular simulations, or estimated using integral equation theo- 
ries. Thermodynamic properties of fluids with pairwise inter- 
actions can be calculated directly from S{k) using exact re- 
lationships from statistical mechanics. Furthermore, many of 
the nontrivial dynamical behaviors of liquids can be predicted 
from a knowledge of S{k) using mode-coupling theory and 
its recent extensions liiSii 

However, despite its considerable practical value, S{k) can- 
not provide a comprehensive description of liquid structure 
mainly because it only contains information about the spatial 
coiTelations between pairs of particles. Higher-order coiTe- 
lation functions, or suitable approximations for them, are re- 
quired to predict structural quantities that depend on the rela- 
tive positions of three or more particles. A well known ex- 
ample of such a quantity is the single-particle free volume 
vi, illustrated in Fig. [2 . It is defined as the cage of acces- 
sible volume that a given particle center could reach from its 
present state if its neighboring particles were fixed in their cur- 
rent configuration. Restated in simple terms, V{ quantifies the 
"breathing room" that a particle has in its local packing envi- 
ronment. The thermodynamic properties of purely athermal 
(i.e., hard-core) fluids can be formally related to the statistical 
geometry of their single-particle free volumes4i^i^ii^ 

The idea that relaxation processes should also be linked 
to free volumes has a long history in studies of the liq- 
uid state,^ '°-" and it has motivated the development of nu- 
merous "free-volume based" models for predicting transport 
coefficients.'^ Unfortunately, a general microscopic frame- 
work has yet to emerge, in part due to the computational and 
experimental difficulties associated with measuring and char- 
acterizing free volumes. However, recent advances in compu- 
tational statistical geometry^ '-' have made it possible to effi- 
ciently calculate such properties from particle configurations 
obtained via either experiments (e.g., confocal microscopy 
of colloidal suspensions'^''') or computer simulations. ^si^iiii^ 
These methods have catalyzed new efforts to quantitatively 
examine the basic ideas underlying the free-volume perspec- 
tive for dynamics and, hence, the prospects for developing a 



successful microscopic theory. 

In this article, we contribute to one aspect of this effort 
by introducing a simple model for predicting the statistical 
geometry of single-particle free volumes in equilibrium flu- 
ids. We use the aforementioned computational tools of sta- 
tistical geometry to test the predictions of the model for (i) 
the hard-sphere (HS) fluid and (ii) a fluid of particles with 
short-ranged, square- well attractions. The former is the stan- 
dard structural reference fluid for simple liquids. The latter 
serves as an elementary model for the behavior of suspen- 
sions of attractive colloids and globular proteins, whose struc- 
tures can be strongly influenced by interparticle attractions. 
As we show, our model suggests a scaling relationship for the 
density-dependent behavior of the HS system. It also predicts 
the manner in which the second virial coefficients of fluids 
with short-range attractions affect their free-volume distribu- 
tions. 



n. GENERAL FRAMEWORK 

We consider a three-dimensional (3D) equilibrium fluid 
comprising N identical spherical particles contained in a 
macroscopic volume V at temperature T. The packing frac- 
tion is 77 = Nira^ /6V, and the particles interact via a short- 
range isotropic pair potential Vij (r) of the generic form 

{00 r < a, 
u{r) a<r<a(l + A), (1) 
r>fT(l + A), 

where A < 1. If one chooses u{r) = — e, then Eq. ([l) 
describes a square-well interaction. Alternatively, the hard- 
sphere potential is recovered if A = 0. We are interested 
in developing a general strategy for predicting how Vij{r), 
1], and T affect the statistical properties of the fluid's single- 
particle free volumes. To simplify notation, we implicitly non- 
dimensionalize all lengths from this point forward by the hard- 
core diameter a of the particles. 

The problem described above cannot be treated exactly, 
and so we make some simplifying approximations that we 
later assess by comparing our predictions to results obtained 
via molecular simulations. Our basic working assumption is 
that the statistical geometry of the free volumes in the 3D 
fluid at T and 77 can be predicted based on a knowledge 
of the exact free-volume distribution /id of the correspond- 
ing one-dimensional (ID) fluid at the same temperature T 
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FIG. 1 : A 2D schematic of the free volume of a tagged particle. Over- 
lapping grey circles represent the exclusion spheres of the neighbor- 
ing particles, (a) The small dashed circle is the tagged particle sur- 
face, and the larger dashed circle is its associated exclusion sphere. 
The cross-hatched region is the tagged particle's free volume, (b) 
An expanded view of the tagged particle's free volume along with its 
approximate representation in our model. 



and the scaled packing fraction ^ = rj/rjMR.T- We choose 
Vmrj = 0.64 to ensure that the close packing limit of the ID 
fluid = 1) maps onto the maximally random jammed state 
of the 3D system.'^ Our specific idea is that a particle's free 
volume in the 3D fluid at T and ry can be modeled as a cuboid 
(see Fig. [3 with length ax{, height ay{, and width azi, where 
X{, j/f , and Zf are independent random variables drawn from 
/id at T and ^. The 3D free-volume distribution f{v{) is thus 



dx dy dz S{vi — a xyz) 
^0 Jo 

X /iD(a;)/iD(y)/iD(^). (2) 



The constant a is a scale factor chosen to ensure that our con- 
struction accurately reproduces the equilibrium free-volume 
distributions and thermodynamic properties of the 3D HS 
fluid (discussed below). Also of interest is the probability 
density (j>{s{) associated with finding a 3D free volume with 
surface area S{. Using our model, 0(sf ) can be expressed as 



0(*f) — / dx dy dzS{s{~2a [xy + xz + yz]) 
Jo Jo Jo 

X flB{x)flD{y)flD{z). (3) 

Eq. (13 and (|3} indicate that, in order to predict f{vi) and 
(j>{s{) using our approach, one only needs to derive the state- 
dependent form of /id from knowledge of Vij (r). However, 
/id is given by 

fmixt) = dyL dyn S{xf - yL - Vr) 
Jo Jo 

X Pg{yL)Pg{yR)dyLdyR, (4) 



where Pg{z) is the probability density associated with finding 
a gap of size z between the surfaces of two neighboring parti- 
cles in the corresponding ID fluid. For particles that interact 
via a pair potential of the form given by Eq. ([0, Gursey2£ has 
shown that Pg{z) can be written 



g-/3n[z+l]g-/3Vij(2+l) 

j^°°g-/3n[.+i]g-/3y,,(s+i)rfs' 



(5) 



where (3 — (fceT")^, and 11 is the ID pressure. The equation 
of state of the ID fluid, i.e., the dependence of H on /3 and ^, 
can be found imphcitly from the foHowing relation)^ 



an 



where 



Jo 



(6) 



(7) 



Having outlined the general strategy of our model, we ex- 
amine some of its predictions for two specific systems in Sec- 
tion III: (i) the equilibrium HS fluid and (ii) an equihbrium 
square-well fluid with short-range attractions. 



III. TESTING THE MODEL 

A. Predictions for the HS fluid 

Here we apply the model of Section II to predict the statisti- 
cal geometry of the free volumes in the equilibrium HS fluid, 
where the pair potential Vij (r) is given by 




(8) 



From Eq. (|8}, (|6j, and (0, it follows that the equation of state 
of the ID HS fluid is given by^i 



/3n/e-i/(i-e). 



(9) 



Substituting Eq. and into Eq. yields the corre- 
sponding gap-size distribution. 



Pg{z) 



1-e 

Then, from Eq. @, /iD(a;f ) is simply 
The first moment of this distribution is 



/iD(a;f) = X[ 



dxf Xffiuixi) 



2;tC/(i-C) 



2(1-0 



(10) 



(11) 



(12) 
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Defining Xf = Xf/{xi), we have 

fiD{xf )dxi = /*D(xf )dXf , 

where 

/rD«)-4.Tf*e-2.;. 



(13) 



(14) 



The main imphcation is that while /id is a function of X{ and 
the packing fraction ^ for the ID HS fluid, fi-[){Xf) can be 
represented by a single curve that is independent of ^. As 
we show below, this property, when used in combination with 
Eq. (13 and (|3} of our model, leads to scaling predictions for 
the density-dependent free-volume and free-surface distribu- 
tions of the 3D HS fluid. 

In particular, making use of Eq. (|2j, (II It . and il2\ . we iden- 
tify that 



(vi) = J^dv vf{v) = a^{xif = {^j^ 
Similarly, combining Eq. (0, dl 1> . and (I12> yields 



(15) 



(16) 



Now, if we normalize the 3D free volume by its first moment. 



V{/{vi), then we have 



f{v[)dv[ = f*{vf)dv{, 



(17) 



The sphericity parameter takes on its minimum possible value 
(Af = 1) for a spherical free volume, while it is larger in mag- 
nitude for less symmetric free volumes. Our model predicts 
that the average value of the sphericity parameter (quantify- 
ing the average shape of the free volumes), given by 



2{xy + xz + yz) 



/•OO f^OG 



X /id(x)/id(2/)/id(^) ^ 1.57, (22) 

is independent of packing fraction for the 3D HS fluid. 

Finally, the equation of state of the 3D HS fluid can 
be formally related to the statistical geometry of its free 
volume s)i^ 



(23) 



Here, P is the pressure of the 3D fluid, and p = N/V is its 
number density. In our model, it is easily shown that 



OO /"OO /"OO 

dx dy dz 
"'0 Jq 



2{xy + xz + yz) 



axyz 



x/id(x)/id(2/)/id(^) = ^^, (24) 



1-e 



and therefore we have 



(25) 



where 



f*{vf) = dx dy dz S{vf — xyz) 
Jo Jo Jo 

X fU^)fmiy)fU^)- (18) 

Similarly, if we introduce the reduced free surface Sf = 
Sf/(sf), then it follows that 



{sf)ds{ = 4i*{sf)dSf, 



(19) 



where 



dx dy dz (5(sf — [xy + xz + yz]/3) 
Jo Jo 

X fm{^)fiB{y)fm{^)- (20) 



Eq. (I18> and (I20> . when viewed together with Eq. il4i . show 
that our model predicts that the scaled free-volume and free- 
surface distributions, f*{vf) and 0*(sf ), are independent of 
packing fraction for the 3D HS fluid. We will return to this 
point later, when we test the predictions of the model via 
molecular simulations. 

To use the model to predict the shapes of the free volumes in 
the 3D HS fluid, we analyze the behavior of the dimensionless 
sphericity parameter Af , defined as^^-^^ 



Af = Sf(67ri/2z;f)-2/3. 



(21) 



B. Simulations of the HS Fluid 

To test the predictions of our free-volume model for the 
3D equilibrium HS fluid, we have performed a series of 
molecular dynamics simulations using a standard event-driven 
algorithm.^"* All runs were carried out in the microcanoni- 
cal ensemble using N — 1000 particles and a periodically- 
replicated cubic simulation cell. Snapshots of the system's 
equilibrium configurations were collected and used to calcu- 
late the geometric properties of single-particle free volumes 
via the exact algorithm presented by Sastry et al.i^ 

Numerical predictions using the free-volume model intro- 
duced here require specifying the value of the geometric scale 
factor a shown in Fig.^ As discussed in Section 11, the goal is 
to choose a in a way that allows the model to provide a reason- 
able overall description of the statistical geometry, and hence 
the thermodynamic properties, of the HS fluid. Throughout 
this work we set a — 0.29, which we obtained from a least- 
squares fit of Eq. ( I15t to the HS simulation data for (vf) versus 
rj over the range 0.375 < i] < 0.525. The quality of this fit is 
illustrated in the inset of Fig.|2t. 

We also compare in Fig. [2^ the free-volume distributions 
f{vf) obtained from simulations with those predicted via 
Eq. (|2j- We observe that the model describes the data rela- 
tively well, with the predictions becoming semi-quantitative 
for the higher packing fractions investigated. This trend with 
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FIG. 2: (a) Free-volume distributions of the 3D HS fluid for pack- 
ing fractions rj = 0.375, 0.4, 0.425, 0.45, 0.475, 0.5, and 0.525. 
Symbols are simulation data and solid curves are the predictions of 
Eq. OJ- The arrow indicates increasing packing fraction rj. Inset: 
Average free volume {v{) versus packing fraction rj. Circles are sim- 
ulation data and the solid line is the prediction of Eq. <15> . (b) Prob- 
ability density associated with observing particles with scaled free 
volume V{ = Vf/{v{). Symbols are simulation data from panel (a) 
and the solid line is the prediction of Eq. <18> . 



packing fraction is expected since the simulated free volumes 
become more compact at higher densities, and, as a result, 
they more closely resemble the simple cuboid shapes that we 
have assumed in our model. 

The same free-volume distributions shown in Fig. [2^ are 
also shown in Fig. I^^, only now scaled in the manner sug- 
gested by Eq. ( I18> . This rescaling yields the probability 
density associated with observing a particle with a particu- 
lar value of Vf = V{/{v{). Interestingly, we observe that, to 
a very good approximation, the scaled distributions obtained 
via simulation for all packing fractions investigated collapse 
onto a single curve predicted by the model. 

As can be seen by Eq. ilOi . the model predicts that the dis- 
tributions of free surfaces for the HS fluid should also collapse 
onto a single curve when scaled in an analogous way. We 
compare in Fig.|3t the scaled free-surface distributions 0* (s|) 
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FIG. 3: (a) Proabability density associated with ob- 
serving particles with scaled free surface area Sf = 
Sf/{s{) in the 3D HS fluid for packing fractions rj = 
0.375, 0.4, 0.425, 0.45, 0.475, 0.49, 0.5, 0.51, and 0.525. 
Symbols are simulation data and solid curves are the predictions 
of Eq. <20> . Inset: Average free surface area (sf) versus packing 
fraction rj. Circles are simulation data and the solid line is the 
prediction of Eq. <16t . (b) Average sphericity (Af ) of free volumes 
in the HS fluid versus packing fraction 77. Circles are simulation 
data, and the solid line is the prediction of Eq. <22> . 



obtained via simulation to the predictions of Eq. ( 120^ Again, 
the collapse of the simulation data is striking, and it is de- 
scribed reasonably well by the free-volume model. The inset 
of Fig.|3^ compares the average free-surface area per particle 
(sf ) obtained from simulations to the prediction of Eq. ( II6I 1. 
As with the predictions for f{vf ), the model provides a more 
accurate description of (sf) at higher packing fractions where 
the free volumes are expected to be more compact. At low 
packing fractions, particles will on average have more nearest 
neighbors which define their free-volume then at high pack- 
ing . Therefore, we expect our cuboid approximation to break 
down at low packing fractions. 

To quantify the compactness of the HS free volumes, we 
compare in Fig. ^jp the average sphericity of the free vol- 
umes (At) calculated from the simulations to that predicted 
by Eq. ( I22t . We find that while the cuboid representation of 
free volumes in ourmodel underestimates (At), it qualitatively 
captures the fact that, on average, the shapes of the free vol- 
umes in the HS fluid are fairly insensitive to changes in pack- 
ing fraction. 

Finally, we test the ability of our free-volume model to pre- 
dict the equation of state of the HS system. In particular, 
we compare Eq. ( I25l l to the well-known Carnahan-Starling2^ 
equation, which provides an accurate description of the HS 
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FIG. 4: Equation of state for the HS fluid as calculated from both 
Eq. <25> and the Carnahan-Starling^^ relationship. Inset: Expanded 
version illustrating that Eq. <25> diverges at t/mrj ~ 0.64, while 
the Camahan-Starling relationship diverges at the unphysically high 
packing fraction of 77 = 1. 



state of the ID square-well fluid can be expressed as 



1 Ae-'^nA^ 



(3U l_e-^/3nA(i_e-/3.)- 



(27) 



Furthermore, by substituting Eq. (|5}, (I26> . and (I27t into 
Eq. 0, one finds that the corresponding ID free-volume dis- 
tribution is given by 

= <^ 2e-'3^(xf - A) - xf + 2A X{&[A,2A), 
' + 2Ae-'3^(l - e-'^') a;f > 2A, 

(28) 



where 



eixf. 



^ g-/3nanf 



1 2 



1 _g-/3nA(i 



(29) 



By substituting Eq. (I28> into Eq. (|2}, one can readily calcu- 
late the free-volume distribution for the 3D square-well fluid 
numerically. 



pressure for packing fractions in the equilibrium fluid range 
< T] < 0.494. As can be seen in Fig. |4] Eq. i25\ shows 
good agreement with the Carnahan-Starling equation for the 
equilibrium HS fluid. Moreover, while the Carnahan-Starling 
equation diverges at the unphysically high packing fraction 
of rj = 1, Eq. j25> predicts that the HS fluid becomes in- 
compressible at the packing fraction of the maximally random 
jammed state (77mr,j — 0.64)i^ 



C. Predictions for the Square- Well Fluid 

To demonstrate that the ideas outlined in Section II can be 
readily extended to other fluid systems, we now use our ap- 
proach to predict the free-volume distributions of a square- 
well fluid with short-range attractions, a basic model sys- 
tem for suspensions of attractive colloids and globular pro- 
teins. Fluids with short-ranged attraction are of particular in- 
terest here because they exhibit free-volume distributions with 
shapes that differ from that of the equilbrium HS fluid at the 
same packing fraction.'** Specifically, interparticle clustering 
induced by the short-range attractions increases the popula- 
tions of both large and small free volumes at the expense of 
mid-sized free volumes. These structural changes play an im- 
portant role in the anomalous dynamical properties exhibited 
by these fluids,"* and the ability to predict such changes con- 
stitutes a sensitive test of our free-volume model. 

The square-well pair potential is given by 

{00 r < 1, 
-e l<r<l + A, (26) 
r > 1 + A. 

Using Eq. i26\ . (|6}, and Q, one can show that the equation of 



D. Simulations of the Square- Well Fluid 

To test the theoretical predictions of our model for the 3D 
square-well fluid, we have performed a series of event-driven 
molecular dynamics simulations^"* using N ~ 1000 parti- 
cles. As with the HS simulations described earlier, a cubic 
simulation cell was employed with periodic boundary condi- 
tions. Equilibrium configurations were stored during these 
runs and later used to calculate^ '^ the geometric properties 
of the single-particle free volumes. For all simulations, the 
range of the square-well attraction was set to A = 0.03. We 
studied the fluid at packing fraction 77 = 0.5 and 0.58. For 
the f] = 0.58 case, a flat distribution of particles sizes with 
mean a and half width S — a/ 10 was used to prevent crys- 
tallization. We examined the behavior of the fluid for various 
values of the strength of the interparticle attraction e, which 
is often quantified using the reduced second virial coefficient 
SI = Here, B2 = (1/2) / dr[l-exp(-/31/,,(|r|))] 

is the second virial coefficient of the fluid of interest, and 
Bf^ = 27r/3 is that of the hard-sphere fluid. For the 3D 
square-well fluid, _B| is given by 

Jo (30) 

= l + (l-e''')[(l + A)3-l]. 

The simulated results for the free-volume distributions are 
displayed in Fig. [5^ and [Sji for 77 ~ 0.5 and 0.58, respec- 
tively. As has been observed for other model fluids with short- 
range attractions,'** decreasing B2 (increasing attractions) in- 
creases the populations of both large and small free volumes 
at the expense of the mid-sized free volumes. These trends 
are captured by the predictions of the free-volume model 
(Fig. 03 and|5ll) using Eq. (|28}, and (|29}. It is evident 
that the model does not match the simulation results as well 
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FIG. 5: Free- volume distributions for the square- well fluid at range 
of attraction A — 0.03, and reduced second virial coefficients ~ 
-0.04, 0.41, 0.60, 0.78, 0.84, 0.91, and 0.94. (a) Simulations at 
packing fraction i] = 0.5 and (b) the free- volume model at described 
in the text at 77 = 0.5. (c) Simulations at packing fraction — 0.58 
and (d) free-volume model at 77 = 0.58. Arrows indicate increasing 



at 7] = 0.58 as at 77 = 0.5. This is most likely due to the 



ID equation of state used in the model being less accurate at 
high packing fractions. Although the free- volume model pre- 
sented here does not reproduce the square-well distributions 
with quanitative accuracy, it clearly captures the important 
physical trends. Moreover, its predictions show surprisingly 
good agreement with the simulation data, given that there are 
no free parameters in the theory that are fit to data for the 
square-well fluid. 



IV. CONCLUSIONS 

In summary, we have introduced a simple model for 
predicting the free-volume distributions of equilibrium flu- 
ids, and we have tested this model using molecular simu- 
lations. The model suggests a new scaling for the density- 
dependencies of the free-volume and free-surface distribu- 
tions of the HS fluid, and these scalings show very good agree- 
ment with simulation results. The free-volume model also 
predicts a reasonably accurate hard-sphere equation of state in 
which the pressure diverges at a packing fraction of 77 = 0.64. 
Finally, the model predicts semi-quantitatively the manner in 
which the attractive strength affects the free-volume distribu- 
tions of a square-well fluid with short-range attractions. These 
considerations suggest that the free-volume model proposed 
here can perhaps be used fruitfully within the free-volume 
based theories for dynamics of fluids to understand how their 
transport coefficients derive from their microscopic interac- 
tions and thermodynamic conditions. 
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